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ABSTRACT 

Context. The observed proportionality between the centripetal contribution of the dynamically insignificant HI gas in the discs of 
£\j , spiral galaxies and the dominant contribution of Dark Matter (DM) - the "Bosma effect" - has been repeatedly mentioned in the 

literature but largely ignored. Since this phenomenology, if statistically significant, tells us something about the relationship between 
the visible baryonic and invisible DM, it is important to re-examine the reality of this effect using formal tests and more modern data. 
Aims. We have re-examined the evidence for the Bosma effect, either by scaling the contribution of the HI gas alone or by using both 
the observed stellar disc and HI gas as proxies. 
Methods. We have calculated Bosma effect models for 17 galaxies in The HI Nearby Galaxy Survey data set. The results are com- 
pared with two models for exotic cold DM: internally consistent cosmological Navarro-Frenk- White (NFW) models with constrained 
compactness parameters, and "Universal Rotation Curve" (URC) models using fully unconstrained Burkert density profiles. 
Results. Fits to spiral galaxy rotation curves computed using just HI scaling are inadequate, despite the clear proportionality seen in 
the outer discs. The poor performance is obviously related to the prominent decrease in the HI surface density in regions of high stellar 
surface density, where HI gas been converted into molecules and stars. The Bosma models that partially correct for this physical effect 
using the stellar discs as additional proxies are statistically nearly as good as the URC models and clearly better than the NFW ones. 
Conclusions. We confirm the correlation between the centripetal effects of DM and that of the interstellar medium of spiral galaxies. 
The edificacy of "maximal disc" models is explained as the natural consequence of "classic" Bosma models which include the stellar 
disc as a proxy in regions of reduced atomic gas. The perception that the Bosma effect could be due to the near-equality of the HI 
surface density and the projected mass density of a cold DM halo is incorrect, both theoretically and empirically. The standard ex- 
planation - that the effect reflects a statistical correlation between the visible and exotic DM - seems highly unlikely, given that the 
geometric forms and hence centripetal signatures of spherical halo and disc components are so different. A literal interpretation of the 
Bosma effect as being due to the presence of significant amounts of disc DM requires a median visible baryon to disc DM ratio of 
about 40%. 

Key words. Galaxies: ISM - Galaxies: kinematics and dynamics - Galaxies: spiral - Galaxies: structure - Cosmology: dark matter 



=5 



O 



o 



- 1—1 

X 



\Q 1. Introduction Given enough interaction between the baryons and CDM ha- 

' los during the formation phase of galaxies, it is claimed that 

\Q One of the pillars supporting the present paradigm of cold this difference can be explained (Navarro, Frenk & White [1996J 

O Dark Matter (hereafter CDM) is the interpretation of the ro- Mashchenko, Couchman & Wadsley 2006; El-Zant, Shlosman 

t-H ; tation curves of spiral galaxies: the latter are thought to re- & Hoffman|200l] Tonini, Lapi & Salucci|2006| Governato et al. 

. quire the presence of radiaUy extended, quasi-spherical a nd very |2010>, but there is still a particular problem for dwarf galaxies, 

^ 1 massive halos of exotic non-baryonic matter (Sanders |20T0T >. since the minimal amount of baryonic mass is unlikely to have 

Nevertheless, there still remain many unanswered questions and been able to affect the dynamically dominant CDM halo, so still 

inconsistencies between theory and observations, particularly a d hoc models fo r the d issipative effects of baryons have to be 

in about the role and effects of the baryons. While the signs of Dark invoked (Oh et al.|2010>. 

Matter (hereafter DM) are particularly evident in massive galax- There are other signs of deep and empirically simple but 
ies with nearly flat rather than Keplerian outer rotation curves theoretically unexpected and largely unexplained connections 
(assuming that the mass-to-light ratio of the disc is not an in- between the visible baryonic and invisible dark components: 
creasing function of radius), the most extreme effect of DM is ac- the Modified Newtonian Dynamics pheno menology (MOND; 
tually seen in dwarf spirals with slowly rising rotation cur ves and Milgro m [1981 ] Sanders & McGaugh 120021 Gentile, Famaey & 
minimal stellar and small gaseous contributions (Moore [T994T deBlok|20T]}; the mass-discrepancy relation (M cGau gh [2004); 
Fits to the rotation curves of both types of galaxies using the ob- the baryonic Tully-Fischer relation (McGaugh 2005 ); and the 
served stellar and gaseous components and various models for universal me an surf ace densities derived both for C DM halos 
the CDM halo have shown that the data are not well-fit using (Gentile et al. |20091 > and baryonic discs (Donato et al. 12009) . 
the theoretically expected Navarro-Frenk- White (Navarro, Frenk The most unusual correlation between the observed proper- 
& White 1997; hereafter NFW) profiles, whose central density ties of disc galaxies and the derived properties of the dark corn- 
profiles are too "cuspy" (deBlok & McGaugh 1997) and simple ponents was identified by Bosma (119781 11981) : the centripetal 
adiabatic contraction of the CDM with the baryons should re- contribution to the global rotation curves by the dynamically in- 
sult in even more cuspy profiles (Sellwood & McGaugh 2005). significant gaseous discs appears to be directly correlated with 



1 



Frederic V. Hessman and Monika Ziebart: The Bosma effect revisited 



that of the remaining DM component - the "Bosma effect". This 
effect has been characterized in various ways by different au- 
thors. After deriving the total disc surface densities implied by 
the rotation curves, Y.d yn (R), and comparing them with the ob- 
served HI or HI+H2 gas surface densities, E gas (R), Bosma ob- 
served "... that as a rule the ratio [of the two] ... is more or less 
constant beyond about one-third of the optical radius, with [HI] 
being the dominant contributor ... in the outer parts" (Bosma 
119811 1. Sancisi d 19991 1 simplified the effect by stating that "... the 
rotation curve in the outer parts can be explained ... by scaling 
up the HI disc by about a factor [of] 10". Hoekstra, van Albada 
& Sancisi (2001, hereafter HvAS) studied the Bosma effect in 
a sample of normal galaxies: they said that "... the total surface 
density of matter needed to explain the observed rotation curves 
... [is] roughly proportional to the surface density of neutral hy- 
drogen" and proceeded to define the proportionality as 

V 2 DM (R) * (|^(*) (D 

HvAS and the others explicitly defined the constant of propor- 
tionality between the centripetal effects of the HI gas and the 
DM as the assumed constant ratio of the DM to HI surface den- 
sities averaged over the same disc, corrected only for the pres- 
ence of Helium. Bosma, however, originally explicitly meant to 
include the total gas surface density, even when he used HI as a 
proxy for the interstellar medium (ISM), so there are potentially 
several different versions of this effect, depending upon what as- 
sumptions are made: 

- what we will call the "simple" Bosma effect or "pure HI- 
scaling" uses only implicitly or explicitly corrected for 
the contribution of He and heavy elements; whereas 

- the "classic" Bosma effect attempts to include the total 
gaseous surface density, either explicitly by using "Lism = 
"Lhi + (again corrected for He and heavy elements), or 
by using the stellar disc as an additional proxy. The latter 
can be done by allowing the effective stellar mass-to-light 
ratio, Tdi SC , to be a free parameter. 

The mean HI Bosma ratio has been measured in almost a 
hundred normal and dwarf spiral galaxies spanning a wide range 
of Hubble types from massive early spirals to late dwarf ir- 
regulars (Bosma 119811 Carignan 119851 Carignan & Beaulieu 
1989; Carignan & P uche [19901 Job in & C arignan [19901 Puche, 
Carignan & Bosma [19901 Swaters [T9991 HvAS; Noordermeer 
2006). A direct correlation between the gravitational contribu- 
tions of the gas and DM is seen in practically all galaxies ana- 
lyzed (Sancisi [1999l Bosma l2004l l. While Swaters ( fT999l l found 
value of 3 -5 (corrected for He) in "late-type" dwarf spiral galax- 
ies, HvAS found a wide range of 5-20 in normal galaxies with a 
median value of about 7 and Noordermeer (2006) a range from 
5-52 and a median value of 17 in early spirals. This behavior 
is supported by a correlation between the value of the Bosma 
factors and the masses of galaxies found by Hoektra et al.: more 
massive galaxies tended to have larger values. With just a few 
exceptions, the rotation curves obtained by scaling up the cen- 
tripetal contribution of the HI gas by some constant factor and 
not including a spherical DM halo resulted in model rotation 
curves which followed the observed curves about as accurately 
as a CDM halo model permits. 

HvAS explicitly set out to test the reality of the Bosma 
effect. After admitting that "...most of the 24 rotation curves can 
be fitted rather well over their full extent", they then proceeded 
to explain it away using various arguments, so that the effect is 



now largely ignored in the literature. A critical review of the 
arguments suggests that their conclusion was premature: 

- "The model curve [of the poorer fits] does not agree with 
the observed rotation curve in the inner region." However, 
they admitted that their treatment of the central bulge, beam 
smearing, and other errors could be the cause of this prob- 
lem. 

- There are "... large wiggles that are not present in the ob- 
served rotation curve." While this is true in a few cases, they 
admit that the wiggles can be created by spiral arm structures 
in which the HI can be converted to H2 and other molecules. 
The HI surface density may, in fact, be a very "noisy" tracer 
of the global gas surface density in galaxies. This is easily il- 
lustrated in the "holes" seen in spatially well-resolved HI im- 
ages of face-on galaxies (e.g. Deul & Den Hartog 1990J. In 
fact, many rotation curves show small-scale structure which 
cannot easily be interpreted as being due to mass-density 
variations using any model. 

- "The model rotation curve drops below the observed rotation 
curve at large radii." The outer scaled-up HI rotation curves 
ended up significantly (k, 10%) below the observed ones in 
about 7-8 of the 24 galaxies, but the deviation always oc- 
curs at the very outer edge of the visible HI disc (93+9%). 
HvAS admitted that there can be several reasons for this ef- 
fect (e.g. large scale flux not seen by the interferometers, 
partial ionization), leading them to conclude that "... the ob- 
served HI surface density may therefore not always trace the 
total gas density in the outer parts." Additionally, the HI ro- 
tation curves are almost always calculated from the observe 
surface densities assuming that the disc is effectively infinite, 
systematically minimizing the rotation curve velocity at the 
edge of the discs (Casertano 1983). 

- "... scaling of HI to represent the dark component only works 
in combination with maximal discs" (Palunas & Williams 
2000). This effect is not surprising, since the models using 
HI alone could simply be correcting for the missing H2 gas 
component most strongly associated with the stars. 

- "... our sample is biased against galaxies with R ou tlhm sub- 
stantially larger than 3." This argument is based on the as- 
sumption that a large disc exponential scalelength is needed 
to insure the flatness of the rotation curves at large radii. 
However, the observed HI distributions in the HvAS sample 
are not particularly exponential: the observe distributions are 
better described as being flat in the inner parts, with nearly 
constant surface densities of 5 - 10 M /pc 2 and having com- 
plex tails (to what extent the total gas surface densities are 
exponential is unknown). They suggested that better obser- 
vations would eventually show the drop in the scaled HI 
curve and hence the failure of the Bosma effect, but Swaters 
( 1999) found no such problem in galaxies where the extents 
of exponential discs are clearly much larger than those of the 
stellar discs. 

- Although HvAS admit that ".. for about two-thirds of the 
galaxies we obtain good fits to the data", one has the im- 
pression that in fact practically all of the galaxies are well fit 
by the Bosma effect. Since the HvAS fits were done by eye 
and were not explicitly compared with standard CDM rota- 
tion curve fits, it is impossible to judge how bad or good a 
"2/3" success rate really is. 

- They finally argue against the reality of the Bosma effect fits, 
saying that "the good fits are somewhat coincidentaF be- 
cause they expect the HI curves to eventually decline very 
rapidly, even though their own data did not, in fact, show 
this. 
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HvAS finally conclude that "... in many cases,, there seems 
to be little or no relation between HI and dark matter" and "... 
it is not possible to conclude here that there is a real coupling 
between HI and dark matter in spiral galaxies." These negative 
conclusions are obviously very conservative: measured on the 
standards of typical rotation curve fits using any mass model, 
the results were, in fact, amazing good. Given the success of the 
phenomenological model - there is no way to modify the shape 
of the centripetal contribution except by global scaling of the 
stars and HI - and a natural explanation for the cases that didn't 
fit quite as well, HvAS could in fact not show the Bosma effect 
is a statistical fluke or selection effect and hence to disprove that 
there was no coupling between the ISM and the gravitational 
effects attributed to DM in spiral galaxies. 

At the time Bosma first described this correlation, placing 
the DM in a disc was a reasonable assumption, but already at 
the time of HvAS, the DM was assumed to reside in a quasi- 
spherical halo. Since the centripetal acceleration V 2 /R at a given 
radius is calculated only indirectly from an integral over the en- 
tire surface mass density for discs and just over the inner vol- 
ume mass density for spherical halos, the connection between 
the surface or volume densities and the rotation curve contribu- 
tions is dramatically different for the two geometries - unlike 
the spherical case, the gravitational effects of matter distributed 
in a disc-like manner effects the motions of test particles both 
within and without the disc, making it much more difficult to 
obtain a particular rotation curve contribution. Ignoring this dis- 
tinction has led to considerable confusion in the literature about 
the role of surface density: e.g. Meuerer & Zheng (2011) in- 
correctly characterize the Bosma effect as showing that "... the 
HI column density in extended HI discs often is nearly a linear 
tracer of the projected DM surface density T The centripetal ef- 
fect of a spherical halo is, in general, not at all the same as that 
of a disc having a surface density equal to the projected density 
of the spherical halo. The normalized ratios of the centripetal 
effects of three different spherical halos and their equivalent pro- 
jected discs are shown in Fig.[T] while there is, of course, some 
similarity in shape, the relative differences for the regions rele- 
vant for the HI discs are about 30%, with gradients in the rela- 
tive effects systematically largest in the inner regions most rele- 
vant to the Bosma effect. The only exception to this rule is that 
the asymptotic projected surface density and V(R) of a singu- 
lar isothermal halo is the same as that of an infinite disc with 
Z(7?) oc R l . However non-singular isothermal CDM halo fits 
yield small but significant core radii (this is one of the main ar- 
guments against "cuspy" NFW profiles; deBlok et al. 2008J, HI 
discs are not infinite and they do not in general show this surface 
density scaling (see Fig.|2j». Thus, the Bosma effect has nothing 
to do with projected CDM density profiles. 

In summary, the usual characterisation of the Bosma effect 
in the literature is wrong on three counts: 
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Fig. 1. A comparison of the centripetal effects of spherical halos 
with compactness c = 10 and the corresponding discs with the 
same vertically projected surface mass-densities. The ratios are 
arbitrarily normalized to the values at R — r core to emphasize the 
similarities in shape. 

with different projected surface densities, and different cen- 
tripetal patterns. 

With this correct description, the Bosma Effekt must be defined 
phenomenologically as 

V 2 DM (R)/Vgas(R) * const (2) 

where the constant can be interpreted as an effective mass-to- 
light ratio contribution above that normally attributable to stan- 
dard ISM components and the origin of the effect is a priori not 
at all obvious. Unfortunately, previous studies have not shown 
how the Bosma effect performs relative to standard CDM mod- 
els. 

If real, the Bosma effect obviously has something to say 
about the relationship between the ISM and DM in disc galax- 
ies. Fortunately, we can now probe the Bosma effect with much 
better means: the quality of the HI rotation curves and surface 
density maps has increased; we can model the stellar popula- 
tions more accurate due our ability to constrain the stellar con- 
tributions in wavelength bands suffering less from interstellar ab- 
sorption; and H2 maps can be constructed from CO observations 
with millimeter interferometers. 

This article is the first in a series looking again at the 
Bosma effect. In this first paper, we compare the performance of 
"simple" (Hl-scaling alone) and "classic" (stellar disc and HI- 
scaling) Bosma effect models with those of standard CDM us- 
ing a set of particularly well-observed galaxies in The HI Nearby 
Galaxy Survey (Walter et al. 120081 hereafter THINGS). 



- the implied physical correlation is between the DM and the 
total ISM, not just its neutral part; the HI has simply been 
used as the most directly observable tracer, despite its ob- 
vious drawbacks, so that the quality of the observed effect 
might be much better (or worse) than previously claimed; 

- the actual comparison is of the centripetal acceleration con- 
tributions of gas and DM at each radius, not the surface 
densities and certainly not the projected surface densities of 
spherical halos; and 

- within the standard CDM paradigm, the implied correlation 
is between two mass components with different geometries, 



2. Rotation curve models 

We fit the data with a standard rotation curve model consisting 
of a stellar disc, a stellar bulge, HI disc, and CDM halo, each 
baryonic component being assigned a mass-to-light ratio T: 

V obs = ^di sc yj isc + ^bulge^luige + "^HlV^ + V^om (3) 

The centripetal effects of the "simple" Bosma effect - pure 
Hl-scaling - are best described by attributing the observed 
atomic hydrogen with an effective mass-to-light ratio 

Tm = 139 + AY ffl = 1.39(1 + f HI ) (4) 
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where the factor of 1.39 corrects the visible gaseous component 
for the presence of He and heavier elements (given the roughness 
of the Bosma model and the uncertainties and systematic errors 
in the HI data, there is no point in worrying about whether a 
correction to primordial or to Solar heavy element abundances is 
more appropriate). Thus, the Bosma contribution to the effective 
mass-to-light ratio, AT///, can be more usefully defined relative 
to the total corresponding mass density of the neutal gaseous 
medium with the scaling factor //// = AT//// 1.39. This is the 
standard definition for the Bosma effect. 

For the "classic" Bosma effect fits, the stellar disc is used as 
a proxy for those inner regions where the density of HI drops 
far below that which one would have extrapolated from outside, 
resulting in an effective stellar disc mass-to-light ratio 



Table 1. Model parameters 



Ydisc = Y* + AT</,- JC = T*(l + fdisc) 



(5) 



The first term is the correction for the true stellar mass-to-light 
ratio and the second, again, is the effect due to using the stars as 
proxies for other - presumedly non-stellar - mass components. 
The reality of AY</ !!t . obviously depends critically on the relia- 
bility of Y*. In the case that Y+ is derived from infrared images 
using detailed stellar population models calibrated by observed 
infrared colours, the true stellar contribution should be fairly re- 
liable. Again, the "classic" Bosma effect measured with the stel- 
lar disc as a proxy is best expressed as a scaling of the total stellar 
mass in the disc, f disc = AY, //sc /Y*. 

In order to quantify the relative success of the Bosma ef- 
fect at explaining rotation curves of disc galaxies with respect 
to normal CDM halo models, we have also fit our data with two 
plausible CDM models: 1) internally consistent theoretical NFW 
halos and 2) the generic Burkert density profiles (Burkert 1995) 
used by "Universal Rotation Curve" assumption (Salucci et al. 
120071 hereafter URC). 

The centripetal contribution of a CDM halo with a NFW den- 
sity profile is 



vl FW (R) 



V 2 

v 20() 

X 



ln(l +c)-^- 



(6) 



where V200 is the rotational velocity at the so-called "virial ra- 
dius", r2oo> where the mean density of the halo reaches a value 
200 times that of the mean cosmic mass-density, the concentra- 
tion factor c = r%x)/r c is the ratio of r^oo to the NFW core scale- 
length r c , and where X = R/r2oo- Normally, cuspy NFW profiles 
are characterized by 2 model parameters; this is the approach 
taken by de Blok et al. (2008 ). However, the N-body calculations 
which resulted in the NFW profile model clearly show that c is 
not an independent parameter but is in fact strongly correlated 
with V200 (or, equivalently, with M200 or r2oo): 



cnfw ~ 7.80 



V200 



100kms~ 



-0.294 



(7) 



(Maccio, Dutton & van den Bosch 120081 their Eqn. 10). 
Including this intrinsic correlation reduces the number of fit pa- 
rameters by one. This approach automatically removes the prob- 
lems with the NFW fits encountered by de Blok et al, who found 
unphysical values of c in 9 of their 20 NFW fits. 

The URC model uses a Burkert density profile, resulting in 



Vmc(R) = 



V. 



200 
X 



ln( 1 + (cX) 2 ) + 2 ln( 1 + cX) - 2 tan~ 1 (cX) 



ln( 1 + c 2 ) + 2 ln( 1 + c) - 2 tan- 1 (c) 



(8) 



In contrast with the NFW model, URC and isothermal halo mod- 
els are usually characterized by the central mass-density p(0) and 



Model 


fdisc 


fm 


CDM parameter(s) 


"Simple" Bosma 


= 


> 


N/A 


"Classic" Bosma 


> 


> 


N/A 


Constrained NFW 


= 


= 


V200 


Burkert (URC) 


= 


= 


^200, c 



the radial scale length of the core, r c . This distinction is unneces- 
sary and even misleading, since the virial radii and concentration 
factors are equally well-defined, making a comparison between 
the models and the fitted parameters more difficult. Thus, we pa- 
rameterize both CDM models with the same quantities, choosing 
the most commonly used variables V200 an d c - 

A summary of the model parameters for each of the four 
models is given in TableQ] Note that each model has very dif- 
ferent constraints: the Bosma models have shapes that are inde- 
pendent of the velocity amplitude (other than the possibilities of 
linear combinations in the "classic" version); the shape of the 
the NFW model is determined by its limiting amplitude; and 
the URC model has no constraints beyond the scaleable generic 
shape. Thus, one naively expects the Bosma models to perform 
worst and the URC to perform best. 



3. The sample 

We use the same subset of THINGS galaxies chosen by de Blok 
et al. (120081 1 for their detailed mass models. Although Gentile, 
Famaey & de Blok (1201 11 1 studied an even smaller subset of the 
THINGS galaxies - those suffering least from non-circular mo- 
tions - we include the full set and adopt the same input data 
(e.g. distances) in order to be able to compare our analyses 
with the unconstrained NFW and isothermal fits performed by 
de Blok et al. The names of the galaxies, their assumed distances, 
NED types, stellar masses, HI masses (multiplied by 1.39 to cor- 
rect for He and heavy elements), fitted exponential radial scale- 
lengths for radii where such are meaningful, and the optical di- 
ameters (B-band R25 from the RC3 catalogue; de Vaucouleurs et 
al. U991l l. converted into kiloparsecs) are given in Table|2] Seven 
of these galaxies were also contained in the HvAS sample, per- 
mitting a direct comparison of the results using better data. 

The stellar contributions for the THINGS rotation curves 
were calculated in deBlok et al. (2008) using 3.6 prni images 
from the Spitzer Space Telescope, thereby reducing the depen- 
dence on ill-constrained extinction corrections. The mass-to- 
light ratios, Y*, were calculated using either a scaled Salpeter 
(Bell & deJong 2001) or Rroupa (1200 11 1 initial mass function 
(IMF): since the scaled Salpeter IMF often resulted in obvious 
overestimates of the disc contribution in the de Blok et al. mod- 
els, we adopt their Kroupa-based mass-to-light ratios as being 
more representative of the true stellar components. The inter- 
ested reader is directed to the detailed discussions in de Blok et 
al. of the handling of colour gradients and bulge components 
when calculating the stellar contributions to the rotation curves 
of individual galaxies. The resulting mass surface density pro- 
files are shown in Fig. [2] as red pluses (when the 3.6 ymi data 
don't extend as far as the HI data, exponential extrapolations 
are shown). The exponential scalelengths of the final disc mass- 
profiles (as opposed to the surface brightness profiles, which are 
described in deBlok et al.) are listed in Table|2] The tabulated 
stellar rotation curves generously provided by the THINGS con- 
sortium already include the effects of the assumed stellar mass- 
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Table 2. The THINGS Sample 



Name 


Type 


D (Mpc) 


Mi, 






Router 


h H i 


for R > R mllel 


lh m 






DD0154 


IB(s)m 


4.3 


0.0019 




0.044 


8.3 


2.3 


2.0 


3.6 


1.9 


1.1 


IC2574 


SAB(s)m 


4.0 


0.074 




0.81 


11.7 


1.2 


7.7 


9.8 


7.7 


2.9 


NGC925 


SAB(s)d 


9.2 


0.72 




0.44 


13.1 




(totally flat) 




14.0 


3.9 


NGC2366 


IB(s)m 


3.4 


0.018 




0.081 


8.2 


1.5 


4.1 


5.5 


4.0 


2.8 


NGC2403 


SAB(s)cd 


3.2 


0.33 


0.029 


0.36 


18.0 


4.9 


7.7 


3.7 


10.2 


1.7 


NGC2841 


SA(r)b 


14.1 


7.6 


1.8 


1.9 


51.7 


15.1 


23.0 


3.4 


16.7 


3.8 


NGC2903 


SAB(rs)bc 


8.9 


1.0 


0.15 


0.89 


29.4 


5.8 


14.0 


5.1 


16.3 


2.5 


NGC2976 


SAc pec 


3.6 


0.13 




0.012 


2.54 


0.7 


1.8 


3.6 


3.1 


0.7 


NGC3031 


SA(s)ab 


3.6 


4.9 


0.91 


0.41 


14.8 


6.3 


12.2 


2.3 


14.1 


2.7 


NGC3198 


SB(rs)c 


13.8 


2.0 


0.20 


2.1 


37.8 


15.0 


10.3 


2.5 


17.1 


3.1 


NGC3521 


SAB(rs)bc 


10.7 


8.7 




1.3 


31.2 


6.1 


17.0 


5.0 


17.1 


4.9 


NGC3621 


SA(s)d 


6.6 


1.4 




1.3 


25.8 


3.7 


17.0 


6.8 


11.8 


2.5 


NGC4736 


(R)SA(r)ab 


4.7 


1.3 


0.36 


0.53 


9.6 


4.8 


2.0 


2.0 


7.7 


1.9 


NGC5055 


SA(rs)bc 


10.1 


8.7 


0.65 


1.4 


48.8 


23.8 


13.0 


2.1 


18.5 


3.9 


NGC6946 


SAB(rs)cd 


5.9 


4.2 


0.27 


0.54 


19.2 


7.4 


9.3 


2.6 


9.9 


3.1 


NGC7331 


SA(s)b 


14.7 


12 


1.2 


1.6 


25.3 


3.4 


20.0 


7.2 


22.4 


4.5 


NGC7793 


SA(s)d 


3.9 


0.19 




1.2 


7.7 


1.7 


5.0 


4.5 


5.3 


1.5 



Notes. Masses in units of 10 10 M G , radii and scalelengths in kpc 



1 /2 

to-light ratios, i.e. are equal to Y+V+iR), so scaling these ve- 
locities immediately yields the Bosma factor fdi SC - 

The T.Hi+He profiles from the THINGS sample (Fig. [2] blue 
exes) along with the stellar profiles (red pluses) - are, in general, 
not exponential, but can be characterized as roughly flat above 
a level of 5 - lOM pc~ 2 , so we have fit an exponential tail to 
the outer regions either below a level of 6 M pc~ 2 (horizontal 
gray lines in Fig.[2j or outside of an obvious central dip at lower 
surface densities ; the corresponding scalelengths are given in 
Tableland displayed in Fig. [2] The ratio of this scalelength and 
the radius of the HI disc varies between 2.0 and 9.8, with a mean 
value of 4.4, so these HI discs extend farther in terms of expo- 
nential scalelengths, than those studied by HvAS. The correla- 
tion between the optical radii, R25, the HI scalelengths, and the 
exponential tail radii, R exp , are shown in the bottom right panel 
of Fig.|2]^e.t/) and hm are typically 90 + 7 and 58 ± 10 percent 
of the optical radii, respectively. Clearly, the break in the HI dis- 
tribution occurs exactly at the optical radius, i.e. where gas has 
been turned into stars. 

Also shown in Fig. [2] are R power-law fits to the HI data: 
the only case where this behavior, correponding to the pro- 
jected density of an isothermal CDM halo, is seen, is perhaps 
in NGC4736. 

The centripetal contributions of the HI gas provided by the 
THINGS consortium were recalculated from the tabulated total 
surface mass densities T.Hi+He without assuming any extrapola- 
tion beyond the last observed radius using the method of Pierens 
& Hure (1990). While this may lead to a slight overestimation 
at large radii (Casertano 1983), it is impossible to be confident 
about how to extrapolate the surface densities even given the at- 
tempts at defining the exponential tails in Fig. [2] so this treatment 
makes the least assumptions about the amount and distribution 
of HI beyond the detection limits. 

We note, as has deBlok et al. (2008, see also Gentile, 
Famaey & deBlok 2011), that the rotation curve error bars do 
not reflect the true statistical errors: they are the result of the 
non-axisymmetries, spiral arms and patchiness of the HI gas dis- 
tribution as well as the difficulties of fitting the HI kinematics 
with a simple model of tilted annulur rings. Thus, the only com- 
parison we can make is to compare our results externally with 
previous analyses and internally by comparing the performance 
of the Bosma effect fits with those assuming CDM. 



4. CDM fits 

We first analysed the THINGS rotation curves by performing the 
constrained NFW and standard URC fits. The resulting parame- 
ters and;f 2 values are shown in Tableland the rotation curves in 
Fig s. [3] H & [5] In addition, we have listed the virial halo masses, 
M200, the fraction of the CDM halo masses contained within (and 
hence directly constrained by) the visible galaxy (defined by the 
last HI rotation curve point), the implied total baryonic mass- 
fractions, and the core radii. Detailed descriptions of the CDM 
fits are given for individual galaxies in Section 7. 

There are clear differences in the median galaxy proper- 
ties between the two different CDM models. For instance, the 
median compactness ratios c are larger (302%), and the me- 
dian CDM masses (41%) and core radii (22%) smaller for 
the Burkert/URC models. Correspondingly, 34% of the median 
CDM mass is contained within the region of the visible galax- 
ies for the Burkert/URC models but only 10% in the constrained 
NFW models. Unsurprisingly, the Burkert/URC fits are clearly 
superior to the constrained NFW fits in almost all cases: the 
median value of the reduced^ 2 values are 3.6 and 1.2 for the 
constrained NFW and Burkert/URC fits, respectively. The most 
extreme failures of the constrained NFW model occur in DM- 
dominated galaxies with very flat rotation curves like NGC2841, 
NGC2903, NGC3031: the correlation between amplitude and 
shape forces the NFW contributions to have too much radial 
gradient to maintain the flat shape. High-mass galaxies like 
NGC7331 have less problems with constrained NFW profiles 
because the massive stellar and fitted DM components can be ar- 
ranged to balance their effects, maintaining the flatness of the ro- 
tation curves. The core radii for the Burkert fits are much smaller 
than those of the NFW ones, resulting in much smaller CDM 
masses. 

While the CDM models generally result in small baryon-to- 
CDM mass ratios, the median value for the constrained NFW 
fits is 4% versus 10% for the URC/Burkert fits. These values 
are consistent with that of the concordance cosmological model 
(20. 1 %) due to the non-inclusion of molecular, ionized, and in- 
falling gas. However, if one prefers the URC/Burkert models, 
one must conclude that these galaxies certainly do not suffer 
much from the "missing baryon problem" (c.f. Bregman 2007) 



5 



Frederic V. Hessman and Monika Ziebart: The Bosma effect revisited 




Fig. 2. Surface densities & exponential fits (straight lines): HI+He (blue exes); stellar disc (red pluses); implied disc DM (black 
diamonds); R l power-law fits to the HI and projected URC profiles (green curves); optical radii R25 (red arrows); radial and surface 
density range for HI exponential scale lengths hui (horizontal and vertial light gray dotted lines). Lower right: comparison of disc 
size with various scalelengths: HI+He (blue exes), stellar disc (red pluses), inner HI cutoff radii (green asterixes). Note that the 
A&A version of this figure has the full resolution. 
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DDQ154 (NFW) 




DDQ154 (Burkert) 



10 12 U 
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NGC2366 (NFW) 
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NGC2403 (NFW) 



NGC2403 (Burkert) 





NGC2403 (Bnsmaj 




NGC7S41 (NFW i 




NGC2841 (Burkert) 




Fig. 3. Constrained NFW (left), Burkert (centre), and "classical" Bosma (right) rotation curve fits. The disc and bulge components 
are plotted as dotted lines, the HI and He contributions as dash-dotted lines, the CDM or Bosma contributions as dashed lines, and 
the total fitted curve as solid lines. Note that the A&A version of this figure has the full resolution. 
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NGC2903(NFW) 



NGC29Q3 (Burkert) 
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NGC303! (Burkert) 



NGCJIvtS (R-jrkcrti 




NGCZ976 (Bosma) 







25 30 35 



NGC3621 (MFW1 




NGC3621 (Burkert) 
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Fig. 4. See caption of Fig. [3] Note that the A&A version of this figure has the full resolution. 
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NGC4736 (NFW) 




NGC.j / .i.ii (Burkert) 




NGC4736 (Bosnia] 




NGCbOSS (Nf-W) 




NGC5055 (Burkert) 





NGC5Q55 (Bosma) 



NGC6946 (NFWj 



NGC6946 (Burkertl 



NGC6946 (Bosma) 






NGC7331 (NFW) 



NGC7331 (Burkert) 



NGC7331 (Bosma) 
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MGC7793 (Bosma) 



NGC7331 (simple Bosmal 





NGC7793 (simple Bosma) 



20 25 BO 




Fig. 5. See caption of Fig. [3] The bottom line shows three "simple" Bosma fits for comparison. Note that the A&A version of this 
figure has the full resolution. 
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and that an apparent dearth of visible baryons can be caused by 
a systematic overestimate of NFW virial masses. 

Interestingly, the CDM fits are not necessarily able to explain 
the last ~10% of the rotation curves: even within the uncon- 
strained Burkert/URC fits, there are many cases in which the fits 
lie significantly either below (DD0154, NGC2841, NGC3621) 
or above (NGC925, NGC3031, NGC7793) the data. Part of this 
problem is undoubtedly due to the difficulty in deriving the mean 
rotation curves for the often corregated or tilted edges (c.f. Oh 
et al. 120081 1. This difficulty was one of HvAS criticisms of the 
Bosma effect - because they did not test the corresponding per- 
formance of CDM models, it was not obvious that this is often a 
generic problem for any model. 

5. The "simple" Bosma effect 

We first checked the Bosma effect by plotting the ratio of the 
centripetal contributions of the DM and gas for each galaxy (the 
"simple" Bosma effect), assuming that the effect is only seen in 
the HI kinematics: 

V DM ^obs ~ T * V di S c ~ r bulgeVl ulge 

~ ~ Hw 2 * fm{ ) (9) 

"gas 1 - Jyy HI 

(see Eqn.[3]). The results, plotted in Fig. [6] show that the inner 
parts of the galaxies have large values of this ratio, i.e. higher 
rotation than can be explained by scaling the centripetal contri- 
bution of HI alone. This is not surprising: the surface density 
distributions of HI in spiral galaxies is known to show flat in- 
ner regions or even severe dips, resulting in small or even neg- 
ative contributions to the centripetal acceleration; and the "sim- 
ple" Bosma model cannot account for the conversion of HI into 
molecules and stars, which preferentially occurs in the centres 
of the discs. 

What is surprising is the relative constancy of the ratio in the 
outer discs. In order to quantify this effect, we fitted the outer 
ratios with a constant, including just enough data to maintain a 
reduced^ 2 < 2.0 (the actual value using realistic errors would 
be lower). The radial regions, fractions of the total disc, and re- 
sulting Bosma factors are listed in Table|4](the factors and radial 
regions are shown as horizontal and vertical lines in Fig. [6]). In 
most cases, the flat regions extend from the optical radius (indi- 
cated by vertical arrows in Fig.|6]l out to the end of the HI disc; a 
formal fit suggests the flat region starts at 92 + 8% of the optical 
radius (the bottom right plot in Fig. [6]), i.e. there is a relatively 
clean "simple" Bosma effect signal once one leaves the region 
dominated by the stars, exactly as originally reported in Bosma 
1981). Unsurprisingly, this radius corresponds exactly with R exp 
(the 8th column in Tableland the bottom right plot in Fig. [2]), 
i.e. where the inner HI surface density profile flattens. This pro- 
portionality exists on the average over 54% of the visible discs, 
with the value ranging from 19% (NGC925) to essentially 100% 
(DD0154). It is important to note that the relative constancy of 
the centripetal ratios between DM and HI exists in the outer discs 
of galaxies independently of the shape of the rotation curve in 
that region: it applies to flat or slightly downward-sloping rota- 
tion curves in massive galaxies as well as to the upward sloping 
rotation curves of dwarf galaxies. 

Given the behavior seen in Fig.|6l it is not surprising that the 
results of formal fits to the entire rotation curves (Table|4| are 
poor, since the "simple" Bosma fit must be made both for the ob- 
viously flat and the obviously non-flat parts of the Vp M /V 2 as ra- 
tio. This behavior is clearly seen for the three representative fits 
shown at the bottom of Fig. [5] The values of fui and reduced^ 2 



are significantly larger than those of the "classic" Bosma models 
described in the next section. Thus, "Hl-scaling" only works in 
conjunction with a proxy for those regions where the HI surface 
density is low due to the presence of stars and molecular gas. 

6. The "classic" Bosma effect 

The results of "classic" Bosma model fits are shown in Figs. [3] 
[4] & [5] and in Table[4] The quality of the fits, as measured by the 
reduced x 1 values relative to those of the URC/Burkert models, 
is somewhat worse, but it is illuminating to qualify the reasons 
on a case-by-case basis; this is done for each galaxy in the next 
section. In any case, the "classic" Bosma fits are better than those 
of the constrained NFW model, even in those cases where the 
latter performs adequately. 

Comparing our values of fm with those of HvAS (Table©, 
one sees that similar values are obtained for NGC2403, 
NGC2841, NGC6946, and NGC7331. The origin of the remain- 
ing discrepancies is two-fold: first, we have exhaustively probed 
the available parameter space while HvAS were only interested 
in finding a solution by hand; and secondly, some of the "classic" 
solutions have values of fdisc which are much higher than those 
which could be normally associated with truly stellar mass-to- 
light ratios, i.e. extreme "maximum disc" solutions. In our case, 
this is not an immediate worry, physically, since we are only us- 
ing the stars as proxies and don't expect the effect to be due to 
an underestimate of the true stellar mass-to-light ratio. 

The values of /^ iJC are generally similar to those seen in 
"minimal disc" fits; the larger values in few galaxies (DD0154, 
NGC2366) are not significant - the stars play a very subordinate 
role in these galaxies. Although the contributions of the stellar 
and HI proxies are different for each galaxy, the median ratio 
of the two is 52%, i.e. the largest contribution is from the "HI- 
scaling" effect. Thus, it is indeed true that the Bosma effect is 
primarily a scaling phenomenology between the ISM and DM 
centripetal contributions. 

7. Results for individual galaxies 

Here is a detailed description of the CDM and Bosma effect fits 
to individual galaxies. 

DD0154 has a rotation curve typical of dwarf galaxies, with 
a rapid rise in the inner disc but no asymptotic flattness. The ob- 
served rotation curve is flatter than the constrained NFW fit for 
radii > 6 kpc and higher than either the Burkert/URC or "clas- 
sic" Bosma fit for the last ~ 10%. Thus, no model is able to fit the 
exact behavior of the outer disc. The "classic" Bosma fit yields 
a very accurate fit to most of the rotation curve, including the 
"dents" in the inner curve which occur at radii of ~l/2, 2, and 
5 kpc. The large nominal value of fdisc is not significant - the 
stars play practically no role in the fit. 

The rotation curve of IC2574 rises continuously over the 
entire disc, like that of DD0154. The constrained NFW fit 
overpredicts the rotation curve at small radii. The mass of the 
Burkert/URC cannot be determined due to the severe degener- 
acy of the model parameters. Both the constrained NFW and 
"classic" Bosma fits fall short of the rotation curve at the outer 
edges. Because of the non-smooth distribution of HI, the Bosma 
fits have too much small-scale structure, resulting in a signifi- 
cant difference in the value of^ 2 relative to the Burkert models. 
Nevertheless, the Bosma fit is still quite good (y v 2 = 1.6). 

The distribution of HI in NGC925 is totally flat and definitely 
cannot be approximated by an exponential. This rotation curve 
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Table 3. CDM Fits 



Constrained NFW (1 -parameter) 





V200 


x? 


c 


loo 


M C DM 


< Rout 




V200 


c 


x, 2 


^200 


McDM 


< Rout 


M vis 
Mcnm 


r c 


DDUid4 


1 a n+0 4 

34.9" 


3.0 


10.6 


50 


1.4 


25% 


ah 


33 


1 n+1 

19 


n a 
0.6 


A £ 

46 


1.1 


33% 


4% 


2.5 


LL,Zj /4 


/i 1 ±1 
41 


0.0 


i n 1 
1U. 1 


39 


2.3 


29% 


3.0 


H 1 +4 


1 n n+0 3 
1U.U 


n 
0.2 


1 n 1 
101 


12 


1 no? 
10% 


HOT 
170 


1 n 
H) 




03 


H 1 
/. 1 


0.9 


90 


0.3 


21% 


1 n 
1U 


1 1 T + l 1 

1 12 


1 n 0+0 1 


1.2 


1 An 
loO 


An 
4 / 


no7 

1 70 


170 


13 


INCjC23od 


QO + 2 

38 


2.3 


10.4 


C A 

54 


1.8 


22% 


5.2 


34 


22 


0.1 


48 


1.3 


34% 


8% 


2.2 


INCjL.24UJ 


139* 1 


2.8 


7.1 


198 


89 


8.7% 


28 


01 1+04 

8 1.1 


no *7+0 4 

28.7 


2.0 


116 


18 


36% 


4% 


4.0 


JN(jL,2o41 


394 


a 1 n 
43.0 


3.2 


303 


21)32 


o.a/o 


1 no 
lOo 


1^/1 + 1 
134" 


42 


2.1 


219 


1 on 
121) 


ZA 07 

34% 


CIO/ 

9% 


3.2 




23o 


42.0 


0.1 


33 / 


43 / 


1.17c 


3d 


1 A/1 + 1 

1U4 


32 


n 
U.O 


1 A O 

14o 


3 / 


32% 


3 70 


1 n 
2.9 


NGC2976 


34 


2.6 


n /i 

9.4 


76 


c 1 

5.1 


2.3% 


8.2 


136!" 


22 


0.5 


1 f\A 

194 


"j 
83 


0.3% 


n on/ 

0.2% 


n n 

9.0 


NGC3031 


278 ±4 


10.2 


5.8 


397 


714 


1.7% 


69 


83* 


46* 1 


3.5 


118 


19 


37% 


33% 


2.6 


NGC3198 


117* 1 


2.3 


7.5 


167 


53 


29% 


22 


94 ±i 


18.9* 04 


1.3 


134 


27 


49% 


16% 


7.1 


NGC3521 


150 ±4 


5.7 


6.9 


215 


113 


16% 


31 


97 ±3 


25* 1 


4.7 


138 


30 


45% 


33% 


5.5 


NGC3621 


126* 1 


0.7 


7.3 


180 


67 


16% 


25 


93* 1 


20.5* 03 


2.3 


133 


27 


37% 


10% 


6.5 


NGC4736 


79 ±2 


3.6 


8.4 


113 


16 


9.1% 


14 


40* 1 


69* 6 


1.4 


58 


2.2 


50% 


102% 


0.8 


NGC5055 


120* 1 


3.9 


7.4 


171 


57 


36% 


23 


133 ±2 


9 9+02 


1.0 


190 


78 


33% 


14% 


19 


NGC6946 


182 ±2 


1.0 


6.5 


260 


200 


5.9% 


40 


104 ±2 


23.0* 05 


1.1 


149 


37 


27% 


13% 


6.5 


NGC7331 


169 ±3 


0.3 


6.7 


241 


160 


10% 


36 


125 ±4 


17 ±i 


0.4 


178 


65 


24% 


23% 


10 


NGC7793 


155 ±3 


4.1 


6.9 


222 


125 


1.8% 


32 


62* 1 


34* 1 


3.1 


89 


8.0 


23% 


18% 


2.6 



URC/Burkert (2-parameter) 



Notes. Velocities in km s , radii in kpc, masses in units of 10 10 M Q , h = 0.7. 
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rises continuously at large scales but shows prominent "bumps" 
at 5 -6 kpc. The latter cannot be explained by any model. Again, 
the mass of the Burkert/URC cannot be determined due to the se- 
vere degeneracy of the model parameters. Both the constrained 
NFW and the "classic" Bosma fits overestimate the curve at 
small radii and underestimate it at larger radii. Correspondingly, 
the Bosma model overpredicts the rotation curve at the disc edge; 
this is exactly the opposite effect than that expected by HvAs. 

NGC2366 has a classic rotation curve that becomes asymp- 
totically flat (though the errors are large for the outer disc). All of 
the models can reproduce the data equally well. As in DD0154, 
the stellar part of the "classic" fit plays no role, so the large value 
of fdisc is not significant. 

The rotation curve in NGC2403 rises rapidly and becomes 
flat/slowing rising for about 3/4 of the visible disc. The con- 
strained NFW fit systematically underpredicts the inner and 



overpredicts the outer halves. The Burkert/URC fit is much bet- 
ter but underpredicts in innermost disc. The "classic" Bosma fit 
does much better in the inner third of the disc, following more of 
the deviations from a simple, smooth curve. None of the models 
can follow the "bumps" in the outer half of the disc. 

The rotation curve of NGC2841 has a maximum around 
10 kpc and very slowly declines out to the edge. Again, the 
constrained NFW fit produces a systematically false slope. The 
Burkert/URC fit does much better, though at the price of a very 
small core radius (10% of the galaxy radius and 1/3 of the HI 
scalelength!) and it also cannot explain the increase in the last 
10% of the disc. The "classic" Bosma model does a much better 
job in the outer disc, roughly follows the shape but shows too 
much structure in the inner 1/2 of the disc, resulting in a slightly 
higher value of^f v 2 = 2.7 compared with the URC fit (2.0). 
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The results for NGC2903 closely follow those of NGC2841: 
the smaller scale wiggles in the "classic" Bosma fit produce the 
big difference in^- v 2 (4.6 versus 0.8) even though Bosma model 
reproduces the basic flat shape fairly well. 

NGC2976 is another dwarf galaxy with a rotation curve 
similar to IC2574. The mass of the Burkert/URC halo cannot 
be determined due to the degeneracy of the model parameters: 
NGC2976 need a roughly linearly increasing centripetal con- 
tribution that can be provided by practically any large value of 
V2oo- All models do a good job of reproducing the rotation curve, 
whereby the Burkert/URC is nominally the best and, again, the 
"classic" Bosma fit shows too much small-scale structure. 

NGC3031 (M81) has a rotation curve very similar to that of 
NGC2841, showing a fairly flat shape accompanied by a glob- 
ally declining level and a big "bump" at ~ 7kpc due to the 
prominent spiral arms. The constrained NFW model predicts an 
extremely flat rotation curve, whereas both the Burkert/URC and 
the "classic" Bosma models can explain the gentle drop to larger 
radii but cannot explain the "bump". 

NGC3198 has a classic rotation curve that rapidly increases 
to a very flat level. The amplitude of the constrained NFW fit 
doesn't permit a small enough core radius to flatten the CDM 
contribution. All of the models have a problem with the bulge, 
but the "classic" Bosma model attempts to compensate by using 
the negative centripetal effect of the HI contribution, resulting in 
an only slightly worse fit that that of the Burkert/URC model. 

NGC3521 is very similar to NGC2841 and NGC2903. The 
constrained NFW model again cannot produce a flat enough con- 
tribution. Both the Burkert/URC and the "classic" Bosma mod- 
els can fit the outer disc but have a difficult time in the inner 10% 
of the disc, resulting in^- v 2 values of 4.7 and 6.1, respectively. 

NGC3621 is very similar to NGC3521. The "classic" Bosma 
fit does a better job in the inner disc but shows too much smaller- 
scale structure in the outer disc, resulting in a higher value ofxv 2 
(4.2 versus 2.2 for Burkert/URC). Both the Burkert/URC and the 
Bosma models underpredict the outer rotation curve. 

The rotation curve of NGC4736 is very similar to that of 
NGC303 1 but has more "bumps" with smaller amplitude. Again, 
the constrained NFW fit is systematically low in the inner and 
high in the outer disc. The Burkert/URC and "classic" Bosma fits 
are essentially equally good, but the former needs a very small 
core radius. 

NGC5055 is very similar to NGC3521: while the CDM 
models do an adequate job, the Bosma models suffer from the 
smaller-scale variations in the HI centripetal contribution due to 
the highly structured HI disc. 

The strong stellar component in NGC6946 permits even the 
constrained NFW model to do a reasonable job at fitting this ro- 
tation curve - all fits have roughly the same value of Xv 1 ( 1 -4, 1 . 1 , 
and 1.3 for the constrained NFW, Burkert, and "classic" Bosma 
models, respectively). Nevertheless, the "kink" at 6 kpc separat- 
ing the inner rise with the outer flat curve is best reproduced by 
the Burkert/URC and "classic" Bosma models. The latter is the 
best at explaining the other "bumps" in the outer half of the disc. 

NGC7331 is very similar to NGC2841, NGC2903, and 
NGC4736. All of the models do a very good job of explaining 
the rotation curve. Only the constrained NFW model is increas- 
ing enough at the outer radius to explain the observed slight in- 
crease at the edge. 

NGC7793 is similar to DD0154, although the rotation curve 
slopes downward in the outer half of the disc, resulting in a poor 
constrained NFW fit. The Burkert/URC model cannot reproduce 
the "dent" at 2 kpc (or, equivalently, the "peaks" at 0.5 and 
4.5 kpc). While the "classic" Bosma model shows too much 



small-scale structure, its ability to follow the medium-scale 
structure results in the best value of^ v 2 (2.4 versus 5.5 and 3.1 
for the CDM models). 

In summary, the generic Burkert/URC model is most suc- 
cessful - as expected. The constrained NFW model is rarely as 
successful as the Burkert/URC or "classic" Bosma models, a sit- 
uation explained by the "cuspiness" of the basic profile as well 
as the strict correlation between the amplitude and the compact- 
ness. The "classic" Bosma model does remarkably well: in most 
cases, the main source of failure of the latter is the presence of 
too much small-scale structure created by the mottled appear- 
ance of the HI disc (spiral arms, etc.). In fact, one has the im- 
pression that an appropriate smoothing of the centripetal con- 
tributions at small scales would result in nearly perfect rotation 
curves in practically all cases. 

8. Discussion 

We can now review the arguments presented by HvAS against 
the reality of the "classic" Bosma effect. 

- "The model curve [of the poorer fits] does not agree with 
the observed rotation curve in the inner region." Our re- 
sults show that the "classic" Bosma model does not have 
this problem any more than the CDM models. On the con- 
trary, the Bosma model is often much more successful in this 
region. Given the fact that we have used the stars as prox- 
ies, this is perhaps not surprising: we have in effect used the 
known efficacy of "maximal disc" models to compensate for 
the negative centripetal contributions of the HI discs, the lat- 
ter produced by the non-exponential HI surface density dis- 
tributions. 

- There are "... large wiggles that are not present in the ob- 
served rotation curve." This effect, when present, is easily 
explained as being due to the inadequacy of using HI as a 
total mass tracer at small spatial scales. On the other hand, 
there are galaxies where the observed "wiggles" in the ro- 
tation curve are more easily explained by the Bosma model, 
where V (R) is more closely tied to the gas density, rather than 
by a CDM model whose small-scale density variations must 
be small (e.g. Springel et al. 2008). 

- "The model rotation curve drops below the observed rotation 
curve at large radii." This effect is obviously present rela- 
tive to the performance of the best CDM model in only 2 of 
the 17 rotation curves (IC2574, NGC2976); the other weaker 
examples show variations in both directions of a magnitude 
corresponding to systematic variations from a smooth rota- 
tion curve and hence not specific to the Bosma model (there 
are also cases where the model lies above the rotation curve). 
The difference between the CDM and Bosma models, while 
present, is not large. If the HI suffers from ionisation and/or 
formation of H2 in the outer disc, this minor effect could 
easily be due only to the ISM chemistry and dynamics in the 
outer, relatively starless disc. 

- "... scaling of HI to represent the dark component only works 
in combination with maximal discs." The flatness of or dips 
in the central HI surface mass-density distributions creates 
centripetal deficits which force the use of a maximal disc 
when fitting "classic" Bosma models. This is obviously due 
to the inadequacy of using the neutral component of the ISM 
alone and is thus expected: the relative success of maximal- 
disc models is a natural consequence of the "classical" 
Bosma effect and not a sign of its failure . 
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- "... our sample is biased against galaxies with R ou t/hni sub- 
stantially larger than 3." This argument does not apply to 
the new data sets used: the mean value of R„ M /hHi for 
the THINGS sample is 4.4, with values between 2 and 10 
present. In any case, this argument by HvAS was originally 
based on the assumption, that the HI is distributed exponen- 
tially, whereas the observed distribution looks more like an 
outer exponential with an inner disc saturated at surface den- 
sities around 5 - lOM pc~ 2 . 

- While HvAS found that ".. for about two-thirds of the galax- 
ies [they] obtain good fits to the data", we find a much higher 
success rate when comparing with the corresponding CDM 
models. Formally, the fraction of fits having values of the 
reduced x 2 less than 2 is 18% for the NFW, 70% for the 
Burkert/URC, and 47% for the "classic" Bosma models. The 
median values of^ v 2 are 3.6, 1.6, and 2.6 for the constrained 
NFW, Burkert/URC, and "classic" Bosma models, respec- 
tively. That the nearly unconstrained URC models perform 
better than the highly constrained Bosma models is, after all, 
not surprising. Given that the systematic errors in the obser- 
vations which produce the sometimes highly structured ro- 
tation curves are not included in the values of x 2 , one can 
conclude that the Bosma model is, in fact, quite adequate. 
Thus, 

- we find that it is not at all true that "the good fits are some- 
what coincidentaF '. 

We conclude that HvAS conservative rejection of the Bosma 
effect was premature and based upon a too subjective judgement 
of the data. While the rotation curve fits are not perfect, they 
hold up very well even against CDM fits with considerable free- 
dom in the amplitude and shape of the DM components - and 
this despite the severe inflexibility of shape present in the "clas- 
sic" Bosma models. Wherever the effect appears to fail, there 
are either usually good reasons why a trivial interpretation of the 
failure is misleading or signs that the problem does not occur for 
the Bosma effect alone, but is a problem for any mass model. 



9. The cause of the Bosma effect? 

Given that the Bosma effect appears to be real, the obvious ques- 
tion is why it occurs. The standard explanation would be that this 
is yet another correlation between the effects of baryons and the 
distribution of CDM. However, the explicit distinction between 
the old (Eqn.Q} and new definitions of the Bosma ratio (Eqn.|2]l 
is not at all subtle: whereas it was originally plausible to think of 
an unknown DM component imbedded within and hence propor- 
tional to the baryonic disc, leading to the definition in Eqn.[TJ the 
simultaneous placement of the DM in a quasi-spherical halo and 
the requirement of direct proportionality of its gravitational ef- 
fects with a dynamically minor component in a disc constitutes a 
difficult problem in an individual galaxy and a seemingly insur- 
mountable problem for a large and dynamically heterogeneous 
sample of galaxies. We have seen that the supposed connection 
between the projected CDM densities and the HI densities men- 
tioned in the literature is a myth based on a misinterpretation 
of the observational basis for the Bosma effect. The observed 
correlation is thus either a tracer of a totally inexplicable inter- 
action between CDM and the baryonic disc or a telling argument 
against the paradigm of CDM in a non-disc halo. 

If one interprects the Bosma effect literally, as being due to 
the presence of disc DM in some form associated with the ISM, 
the amount of additional mass is significant but not totally un- 



reasonable. The total implied disc DM masses, Mjhm, and the 
corresponding visible-to-disc DM ratios 

Myis _ M disc + M bu i ge + M HI+He 

MdDM fdiscMdisc + fm^HI+He 

in Table|4] suggest that the visible mass makes up 11-63% of 
the total, with a mean of 25%. The derived distributions of the 
Bosma component are also shown in Fig. [2] (black diamonds): 
one has the striking impression that they are exponential with 
large scalelengths more often related to the outer HI disc than to 
the stellar disc. The "classic" Bosma effect appears to conspire 
so that the dent in an exponential gaseous disc visible in the HI is 
filled up by a scaled version of the stellar disc, revealing a glob- 
ally exponential distribution of disc DM. Also shown in Fig. [2] 
are the projected surface densities of the Burkert/URC fits: with 
the exception of NGC2841, the disc DM - if it exists - is not 
distributed like the projected CDM halos. 

A candidate for disc DM that would naturally explain the 
Bosma effect is a nearly invisible ISM component consisting 
of very dense and cold H2 cloudlets (Pfenniger, Combes & 
Marti net [19941 Pfenniger & Combes [19941 Gerhard & Silk 
1996) whose dynamical effects within the galaxies could be 
very different from those usually assumed for the visible mat- 
ter (Revaz et al. |2009t . Indeed, there are many reasons based on 
the formation and dissociation of H2 to believe that the visible 
HI structures in galaxies belie only a fraction of the true gaseous 
content, fully independently of any questions concerning their 
centripetal signatures (Allen 2004). On the other hand, it is gen- 
erally argued that it is impossible for the discs of spiral galaxies 
to contain enough matter to explain the rotation curves for a va- 
riety of reasons, making it easy to dismiss this interpretation of 
the Bosma effect on largely theoretical grounds. This is not the 
adequate place to address each of these arguments in detail - 
this will be the subject of a future paper in this series - but it 
is clear that they are based upon sets of assumptions or inter- 
pretations that may sound reasonable but are not guaranteed to 
be correct, at least as far as the necessity of exotic cold DM is 
concerned (warm or hot DM cannot be adequately probed at the 
small scales used here: Sanders 120071 Gentile, Zhao & Famaey 
2008). In any case, the mass needed to explain the Bosma effect 
using disc DM is still much lower than the amount of "missing 
baryons" within the concordance ACDM model (c.f. Bregman 
120071) . so the question is not one of having enough baryons but 
whether one can place them within stable galaxy discs in a form 
that is not easily observed. 

If the Bosma effect is telling us that the baryons have sim- 
ply not been given enough weight, then the only other alter- 
natives are a change in the properties of gravity like MOND 
(see Gentile, Famaey & de Blok 1201 II for an analysis of the 
same data) and Conformal Gravity (Mannheim & O'Brien BOl II) 
or the presence of non-linear effects in self-gravitating systems 
within General Relativity (Cooperstock & Tieu 2005 ; Carrick & 
Cooperstock l201 ll but see Fuchs & Phleps 2006), possibly re- 
sulting in non-Newtonian behavior even in low-gravity systems. 
Which of these different explanations one finds objectionable, 
plausible or appealing is, of course, a matter of perspective. 

The Bosma effect is a significant addition to the other known 
but poorly understood correlations between the baryonic com- 
ponents of spiral galaxies and their DM components. As such, 
it deserves much more study. In following papers, we will make 
a better correction for the amount of matter in the ISM in order 
to further constrain the origin of this effect and present a com- 
parison of the Bosma effect with other baryon-DM correlations 
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Fig. 6. Ratio of the centripetal contributions of the DM and HI gas as a function of radius (red horizontal line), corresponding to 
1 +/hi in the "simple" Bosma model. The fitted radial region was selected between the vertical blue lines. Bottom right: optical radii 
(.R25) versus the starting radius where the ratio is statistically flat. Note that the A&A version of this figure has the full resolution. 
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